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ABSTRACT 


The land surfaces of hot-humid tropical urban areas are exposed to significant levels of solar 
radiation. Increased heat gain adds to different land surface temperature profiles in cities, 
resulting in different thermal discomfort thresholds. Using multi-temporal (1986, 2001, and 2017) 
landsat data, this study examined the impact of land use change on urban temperature profiles in 
Umuahia, Nigeria. The findings revealed that over time, built-up regions grow in surface area and 
temperature at the expense of other land use. The transfer matrix, showed that approximately 59.88 
percent of vegetation and 8.23 percent of bareland were respectively changed into built up during 
the course of 31 years. The highest annual mean temperature in built-up regions was 21.50°C in 
1986, 22.20°C in 2001, and 26.01°C in 2017. Transect profiles across the landuses reveals that 
surface Temperature rises slowly around water/vegetation and quickly over built-up and bare land 
area. The study observed drastic changes in land cover with a corresponding increase in surface 
temperature for the period between 1986 and 2017 with consistent decrease in water bodies and 
bare land in the study area. Overall, the spatio-temporal distribution of surface temperature in 
densely built up areas was higher than the adjacent rural surroundings, which is evidence of Urban 
Heat Island. The impact of landuse change on urban surface temperature profiles could provide 
detailed data to planners and decision makers in evaluating thermal comfort levels and other risk 
considerations in the study area. 
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1.0. Introduction 


Different land uses in urban areas result in a variety of natural and man-made surfaces which alters 
their surface structure, energy balance, and micro-atmospheric conditions. Urban development in the 
humid topics is marked by the movement of people from rural areas to cities (Badmos et al., 2020; 
Tacoli, 2001). The footprints of urbanization naturally decrease vegetation cover and modify the land 
surface to the point where local and regional temperatures are influenced (Pielke, 2005; Sewall et al., 
2000). Temperature profiles of cities in the humid tropics are modified by anthropogenic factors 
resulting from many physical and economic activities. These changes in temperature can adversely 
impact our health, occupation, infrastructure and thermal comfort (Edward et al., 2021). 


Studies such as Arnfield (2003) and Mukherjee and Singh (2020) have found considerable thermal 
differences in urban areas, particularly during the day, as well as a strong association between land 
use changes and temperature distribution. These shifts are linked to changes in urban land use and 
land cover patterns including physical and morphological factors. 


It's not unexpected, then, that human activity in the process of urbanization decreases vegetation cover 
and alters land surface to the point that regional and global climates are altered. As a result of these 
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changes, surface temperatures rise several degrees higher than temperatures in adjacent rural areas, 
resulting in the urban-heat-island (UHD effect. UHIs are largely controlled by background climate, 
such as cloud cover (Zhao et al., 2014) and city size (Oke, 1973). Unfortunately real-time records of 
surface temperature changes due to land use change are missing in most cities of Nigeria. In 1960, the 
population of Umuahia urban was only 23,000 people; however, by 2021, it had grown to 816,983 
people, a 5.58 percent increase. Furthermore, the recent relocation of the main commercial complex 
from the city center and the construction of various administrative buildings illustrate Umuahia's 
urbanization patterns. The spatial information relating to these land-cover change-temperature 
relationships is lacking in the study area. 


The dearth of continuous point surface temperature data in Nigeria's south-eastern region is owing to 
a poor network of meteorological stations. Individual measurements obtained from these point data 
sets, are not a fair depiction of the spatial and temporal profiles of temperature in reality. Remotely 
sense Land Surface Temperature provides pertinent information about surface physical properties, 
energy balance (Hu et al., 2020; Mallick et al., 2020; Tran et al., 2006) and heat diffusion in 
landscapes. Thermal remote sensing allows for the measurement, analysis, and modelling of land 
surface temperature, which is a key aspect in a variety of environmental processes such as global 
warming and UHI (Abdulla - Al Kafy et al., 2021). Variations in temperature across different land 
uses have the potential to affect the general landscape and human activities. 


Diverse landscape morphology and human footprints in cities create complex micro-climatic 
conditions and surface temperature profiles in the humid tropics. Due to minimal yearly diurnal 
temperature changes in the tropics, these conditions exist all year. The minimal diurnal temperature 
variability makes implementing climate-adaptable measures into human activities difficult. 
Understanding the relationship between the urban landscape and its impact on the microclimate, on 
the other hand, could be valuable in defining mitigation elements that could greatly contribute to 
achieving a comfortable thermal environment (Sharifah ef al., 2012). The aim of this study is to 
provide understanding on the influence of land use change and its impact on surface temperature 
profiles in a hot humid urban area. This is conducted through the use of remote sensing and field 
investigation in Umuahia, South-Eastern Nigeria. 


2.0. Methodology 


2.1. The study area 


Umuahia is located in the South-Eastern part of Nigeria along Latitudes 5° 25’0” to 5° 4330” N and 
Longitude 7° 22’0” to 7° 36’0” E. The land area is about 385.02 km2 with a population of about one 
million people (Figure 1). The climate is humid tropical rain forest, with a dry season (November- 
March) and a rainy season (April-October). Umuahia's climate is classified as typical hot and humid 
tropics, with warm, sunny days, consistent temperatures, and high humidity. The average annual 
rainfall is approximately 2500 mm. 


Umuahia is drained by the Imo River on the west axis, the Kwa-Ibo River on the south axis, and 
Enyong Creek on the east axis. River Eme empties into the Imo River, while the Ofenyi River 
meanders through the Ibeku countryside before emptying into the Enyong River. 


2.2. Data 


Freely available Landsat images for 1986 (TM), 2001 (ETM+) and 2017 (OLI) were downloaded 
from https://earthexplorer.usgs.gov and processed to extract land cover and temperature matrices for 
the study area (Table 1). Geometric correction of the images was carried out using 241 ground control 
points (GCP) collected through Global Positioning System (GPS), and a 1967 topographic map 
(1:50,000) covering the study area. The selection criteria for the images involved three criteria, 
namely: (1) The time difference between the images was at least five years; (2) the images had similar 
air temperature when acquired; and (3) the datasets were cloud free (Deilami et al., 2016). 
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Figure 1: Location of the study area 
Table 1: Sources of Satellite Data 

Satellite Sensor Path/Row | Date Data type & Bands | Source of Data 
LANDSAT-S5 TM. 056/188 12-19-1986 | Digital (2-,5) USGS 
LANDSAT-7 ETM+ 056/188 12-21-2001 Digital (1-7) USGS 
LANDSAT-8 OLI/TIRS 056/188 12-26-2017 Digital (1-11) USGS 


























2.3. Data analysis 


2.3.1. Image processing and extraction of surface parameters 

To retrieve surface temperature and land cover types from the image bands, the images were 
processed to derive the Digital Number (DN) and the Top-Of-Atmosphere (TOA) values. The 
radiance values were then used to derive the underlying land cover patterns and LST. The Landsat 
TM/ETM, digital numbers (DNs) were converted into TOA radiance and reflectance using Equations 
1 and 2. 


_ LMAX, — LMIN) 


mn x (Qeal = Qcaimin) + LMIN) 


= (1) 
Qcalmax a Qcalmin 
_ (wx, xd?) (2) 
4 (ESUN, x sin(®) 
L, is the spectral radiance at the sensor aperture in W/(m° sr um); Qa is a quantized calibrated digital 
number; Qeaimax 20d Qeaimin ate the maximum and minimum quantized calibrated pixel values derived 
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from the metadata file of the images. LMAX, and LMIN, are the spectral radiance scales of Qeatma, and 
Qeatmins a represents the distance between the Earth and the Sun in astronomical units; ESUN, equals 
solar irradiance and @ is the sun’s elevation in degrees. 


For Landsat OLI Data 2017, DNs of OLI image for reflective (2, 3, 4, 5, 6, and 7) and thermal (band 
10) bands were converted into TOA radiance and surface reflectance based on Equations 3 and 4. 


LA = MLQcal + Al (3) 
_ PA 4 
7 sin(Ocp) 4) 


LA is the TOA spectral radiance, MLBand-specific multiplicative rescaling factor. Qi is quantized 
and calibrated standard product pixel values, and Al is the band-specific additive rescaling factor from 
the metadata. PA is the corrected TOA reflectance for the sun’s angle; and refers to the sun’s elevation 
angle in degrees, derived from the metadata file associated with the OLI image. 


2.3.2. Land cover classification and accuracy assessment 

Four land cover classes were generated, namely; built-up (including paved roads, residential and 
factory buildings), bare soil, water bodies, and vegetation were generated using maximum likelihood 
classification. For each image, 100 samples were selected, and 25 samples were given to each land 
cover class. The data containing all classes was collected from homogenous areas during field survey. 
Hence, in selecting the training strategy, many factors were considered, including the number of 
pixels used, size of the data, the impact of spatial autocorrelation, and differences in the image. 

In this study, 80 random points were selected for accuracy assessment using a 1967 Top sheet of 
Umuahia (1:50:000), GPS acquired ground control points and Google maps to check the accuracy of 
the classified maps. The overall accuracy of the land cover maps for 1986, 2001, and 2017 was 98%, 
94 %, and 98%, respectively (Table 2). 


2.3.3. Derivation of Brightness Temperature (BT) from the satellite images 

The Brightness Temperature from TM and OLI were estimated using thermal bands 6 and 10 
respectively, based on (Avdan and Jovanovska, 2016). Accordingly, the TOA radiance of the thermal 
was converted into at-sensor brightness temperatures (BT) using the thermal constants given in the 
metadata file (Table 3) and Equation 5. 


Br = in +1) (5) 


Br is the at-sensor brightness temperature in Kelvin; L; equals the TOA radiance. K; and K> are pre- 
launch calibration constants. 


2.3.4. Retrieving land surface emissivity and temperature 

Land Surface Emissivity (LSE) is a proportionality factor that scales the black body radiance (Plank’s 
law) to measure emitted radiance and it is the ability to transmit thermal energy across the surface into 
the atmosphere (Avdan and Jovanovska, 2016). At the pixel scale, natural surfaces are heterogeneous 
in terms of variation in LSE. In addition, the LSE is largely dependent on the surface roughness and 
the nature of vegetation cover (Mallick et al., 2008). These can be corrected according to the nature of 
land cover by deriving corresponding emissivity values from the Normalized Differences Vegetation 
Index (NDVI) values for each pixel. The Emissivity was derived using Equation 6. 


eX = evAPv + esA(1 — Pv) + CA (6) 


€, and €, are the vegetation and soil emissivity respectively, and C is the surface roughness taken as a 
constant value of 0.005 (Jiménez-Mufioz and Sobrino, 2003). 


NDVI —NDVI¢ \2 
v= ( =.) a) 


NDVly — NDVI, 
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P,, is the proportional vegetation. NDVI is the minimum and NDVIv equals the maximum. The 


generated brightness temperature values were then converted to emissivity-corrected LST using 
Equation 8 (Artis and Carnahan, 1982). 


Br 
LST = ——__—__._ 8 

1+ (AoTg / (hc)) Ine (8) 
4 is the average wavelength, is the Boltzmann constant (1.38x10~ J/K), A is Plank’s constant 


(6.626x10™* Js), c is the velocity of light in a vacuum (2.998x108 m/s) and is the emissivity ranges 
between 0.991 and 0.973). Finally, the derived LSTs were converted to degrees Celsius from Kelvin. 


2.3.5. Change detection 


Basic image subtraction and cross tabulations were used to detect changes and transformations in the 
land cover and temperature pixels. Land conversion matrices were also generated in order to 
qualitatively and quantitatively report on the modes of change in the temperature and land use classes. 
The software used in this study includes: Arc GIS 10.4, ERDAS imagines 9.1 and ILWIS 3.3. In 


addition, for the purpose of conducting analysis and creating charts and graphs, Microsoft Excel 2013 
was utilized. 


Table 2: Accuracy assessment for the land use types 
































































































































3.0. Results and Discussion 


3.1. The spatial pattern of land cover change and conversion 


Land cover maps generated for 1986, 2001, and 2017 are presented in Figure 2 (a-c). As at 1986, 
vegetation had the highest spatial coverage followed by bare land, built up and water bodies. 
However, by 2017, built up had increased by 21% at the expense of other land (Table 4). 
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Land Water Built-up Vegetation Bare land Row Total User 
Cover Accuracy (%) 
B S 5 B S 5 B S q 2 S 5 2 S 5 2 S 5 
a S S ro S S a S S x S S a S S a S S 
onl N nN baal N N = nN N _ N N se) N N coal N N 
Water 20 18 19 0 0 0 0 0 0 0 1 0 20 19 19 100 95 100 
Built-up 1 1 20 20 20 0 0 0 0 1 0 20 22 21 100 91 95 
Vegetation 0 1 0 0 0 0 19 20 20 0 1 0 19 22 20 100 91 100 
Bare land 0 0 0 0 0 0 1 0 0 20 17 20 21 17 20 95 100 100 
Column 20 20 20 20 20 20 20 20 20 20 20 20 80 80 80 
Total 
Producer 100 90 95 100 100 100 95 100 100 100 85 100 
Accuracy 
(%) 
Year 1986 | 2001 | 2017 
Overall 98 94 98 
Accuracy 
(%) 
Table 3: Metadata thermal values of the images 
Variable Description Value for | Value for Value for 
landsat 5 | landsat 7 landsat 8 
™ ETM+ OLI 
K, Thermal constants, 607.76 666.09 774.8853 
K, 1260.56 1282.71 1321.0789 
Lmax Maximum and Minimum values of | 15.303 12.650 22.00180 
Lmin Radiance 1.238 3.200 0.10033 
Qcalmax Maximum and Minimum values of | 255.0 255.0 65535 
Qcalmin Quantize Calibration 1.0 1.0 1 
Band Thermal bands Band 6 Band 6 Band 10 
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Figure 2: Changes in land cover for (a) 1986, (b) 2001 and (c) 2017 


Table 4: Changes in land cover categories 
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Land Cover Total Area Coverage in Hectares % Change 

1986 (%) 2001 (%) 2017 (%) 1986-2001 | 2001-2017 | 1986-2017 
Bare land 4571.93 (11.9) 2229.06 (5.8) 484.76 (1.3) -6.1 -4.5 -10.6 
Built-up 2900.66 (7.5) 5363.84 (13.9) | 10978.44 (28.5) | 64 14.6 21 
Vegetation 30905.68 (80.3) | 30805.80 (80.1) | 26994.20 (70.2) | -0.2 -9.9 -10.1 
Water 100.49 (0.3) 80.06 (0.2) 21.37 (0.1) -0.1 -0.1 -0.2 











The transfer matrix (Table 5) reveals different types of transfer probabilities while quantitatively 
demonstrating the land cover transfer process. The rows in the table signify the land cover status and 
transfer the situation in the preliminary state (1986) of land cover change, while the columns of the 
table represent the land cover status and transfer the situation in the final state. 


Notably, a huge portion of vegetation and bare land were converted to built-up. In addition, water 
bodies and bare land decreased considerably between 1986 and 2017 by almost 10.8%. The massive 
conversion of other land cover types (Figure 3) reported in this study is in agreement with similar 
results of (Musakwa and Van Niekerk, 2013; Adewole et al., 2020). Built-up areas in cities in 
developing countries are increasing more than other land cover types. 
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Figure 3: Detected land use conversion (1986-2017) 


Hence, the conversion of vegetation to built-up in Umuahia could be attributed to changes in the 
physical landscape, namely: the increase in administrative structures such as the new government 
house, and Okpara square market, and the recent relocation of certain trade centres within the study 
area. 


Table 5: Land cover transformation matrix 












































Land Cover 2017 
3S | Class Built-up Bare land | Vegetation | Water | Grand Total (ha) 
S| Built-up 3453.76 173.23 732.94 0.96 4360.89 
3 | Bare land 898.66 46.70 479.30 - 1424.67 
6 |_Vegetation 6540.46 275.08 25800.44 2:95 32621.52 
2 | Water 29.81 0.17 40.00 1.71 71.69 
S| Grand Total (ha) | 10922.69 | 495.19 27052.68 8.22 38478.77 








In order to visualize temperature trends on each of the land cover types, 3 transect were generated 
(Figure 4a-c). Transects A (the northernmost, which follows a west to east line starting at Umu Agu 
community), B (a mid-city transect that runs from the edge of the Imo River, moving through Okahia 
community and Amuzukwu and ending at Avonkwu) and C (from the boundary between the Imo 
River and Umu Nwawa, then runs through Avodim to Okwu community). All of the transect profiles 
demonstrate a pattern of lower temperatures near water and vegetation surfaces, increasing sharply 
along built-up and barren land sections. The temperature variations found in each of the land covers 
studied are similar to those found in other studies (Zi-Qi et al., 2017; Magidi and Ahmed, 2020). It 
further reveals that impervious surfaces had a greater surface temperature than non-impervious /water 
surfaces. 
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Figure 4(a-c): Transects A ( the northernmost which follows a west to east line starting at Umu Agu 
community), B (A mid-city transect that runs from the edge of Imo River, moving through Okahia 
community and Amuzukwu and ended at Avonkwu) and C (from the boundary between Imo River 

and Umu Nwawa then runs through Avodim to Okwu community) 


3.2. LST changes on land cover 


Table 6 shows that the average surface temperature increased by 0.43 °C between 1986 and 2001, and 
by 1.92 °C between 2001 and 2017. The spatio-temporal distribution of surface temperature shows 
that built-up areas recorded the highest mean temperature followed by bare land, water and 
vegetation. Vegetation showed a considerably lower radiant temperature in the first two periods 
because it probably reduced the amount of heat stored in the soil and surface through transpiration. 
The surface temperature of the water was generally low when compared to other classes. However, 
the increase rate was highest throughout the period (Table 7). The maximum decadal change in 
temperature surface temperature for the different land cover was recorded for water (1986-2001) and 
Built-up/bare land (2001-2017). Correspondingly, built-up accounted for the maximum temperature 
increase for the 31 year (1986-2017) period. 


440 Ike et al., 2021 


Nigerian Journal of Environmental Sciences and Technology (NIJEST) Vol 5, No. 2 October 2021, pp 433 - 443 


Table 6: Mean land surface temperature 










































































Land 1986 Land Surface Temperature (°C) 2001 Land Surface Temperature (°C) 2017 Land Surface Temperature (°C) 
cover 
classes 
Area (ha) Min Max Mean Area (ha) Min Max Mean Area (ha) Min Max Mean 
Bare land | 4571.93 19.12. | 23.61 [ 21.37 | 2229.06 19.44 | 23.55 | 21.50 | 484.76 22.61 | 28.00 | 25.31 
Built-up 2900.66 19.10 | 23.90 | 21.50 | 5363.84 20.84 | 23.55 | 22.20 | 1097844 | 23.30 | 28.72 | 26.01 
Vegetation | 30905.68 18.08 20.40 19.24 30805.80 15.94 22.96 19.45 26994.20 20.78 25.53 23.15 
Water 100.49 19.28 | 20.04 | 19.66 | 80.06 19.54 | 22.81 | 21.17 | 21.37 21.46 | 22.81 | 22.13 
Table 7: Nature of surface temperature increase in different land cover classes 
Class Increase | Yearly | Decadal | Increase | Yearly | Decadal | Increase | Yearly | Decadal 
in °C increase | increase | in °C increase | increase | in °C increase | increase 
from from from 
1986- 2001- 1986- 
2001 2017 2017 
Bare land | 0.13 0.01 0.09 3.81 0.24 2.38 3.94 0.13 1.27 
Built-up 0.7 0.05 0.47 3.81 0.24 2.38 4.51 0.15 1.45 
Vegetation | 0.21 0.01 0.14 3.7 0.23 2.31 3.91 0.13 1.26 
Water 1.51 0.10 1.01 0.96 0.06 0.60 2.47 0.08 0.80 






































The land cover maps revealed that the surface temperature had increased at a high rate in those areas 
where other land cover classes were converted to built-up (Figure 5a-c). It was also observed that 
temperatures continued to increase overtime across the study area in response to the increase in built- 
up areas and bare surfaces. The continuous human-induced transformation of the existing land cover 
to built-up impervious surfaces increases the albedo of the surface and hence increases land surface 
temperature around the city centre of Umuahia. Similar results have been reported elsewhere 
(Kottmeier et al., 2007; Sahana et al., 2019) as natural surfaces have been found to exhibit lower 
temperatures than built-up. 


























Figure 5: Land Surface Temperature Map of (a) 1986, (b) 2001 and (c) 2017 


4.0. Conclusions 





The study observed drastic changes in land cover with a corresponding increase in surface 
temperature for the period between 1986 and 2017 in Umuahia, South-Eastern Nigeria. Land 
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conversion matrices also revealed a systematic conversion of vegetation into built-up and bare land. 
Similarly, there was a consistent decrease in water bodies and bare land for the period under review. 
Overall, Built up recorded the highest temperature records followed by bare land, vegetation, and 
water bodies. Overall, the spatiotemporal distribution of surface temperature in Umuahia's built-up 
areas was higher than in other land uses, indicating the possibility of an Urban Health Island. Further 
studies on the Intensity of UHI (e.g. Bassett et al., 2021) are required to substantiate ours claims. 


The impact of landuse change on urban surface temperature profiles could provide detailed data to 
planners and decision makers in evaluating thermal comfort levels and other risk considerations. 
Despite the fact that a better understanding of landscape influences on the thermal environment may 
aid practitioners in establishing mitigation measures, obtaining outdoor thermal comfort in a hot- 
humid climate is not an easy task. As a result, research into hot-humid temperature profiles, 
particularly the ability of natural vegetation and hydrological factors to regulate local climate, should 
be expanded and enhanced in order to get a better understanding of the individual climatic needs. 
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